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We study the influence of the pulse duration on high harmonic generation (HHG) with exploring a 
wide laser-parameter region theoretically. Previous studies have showed that for high laser intensities 
near to the saturation ionization intensity, the HHG inversion efficiency is higher for shorter pulses 
since the ground-state depletion is weaker in the latter. Surprisingly, our simulations show this high 
efficiency also appears even for a moderate laser intensity at which the ionization is not strong. A 
classical effect relating to shorter travel distances of the rescattering electron in shorter pulses, is 
found to contribute importantly to this high efficiency. The effect can be amplified significantly as 
a two-color laser field is used, suggesting an effective approach for increasing the HHG yield. 

PACS numbers: 42.65.Ky, 32.80.Rm 


I. INTRODUCTION 

Because of the promising application as the attosecond 
light source [Il-Q, high harmonic generation (HHG) has 
attracted great interests in recent years M- According 
to the well-known three-step model Q, the maximal en¬ 
ergy of the harmonic (the cutoff energy) in the HHG is 
Ip + 3.17C/p. Here, Up = Eq/{Au)q) is the ponderomotive 
energy with Eq and wq being the laser amplitude and 
frequency. Ip is the ionization potential of the ground 
state. To obtain shorter attosecond pulses, higher cutoff 
energy is expected. To increase the HHG cutoff, one can 
increase the laser intensity or the wavelength. However, 
there are limitations for both manipulations. First, high 
laser intensities can induce the important depletion of 
the ground state which decreases the HHG yield [1, 0 . 
Secondly, due to the diffusion effect, the HHG yield also 
decreases very fast as the laser wavelength increases [13- 
[T3| . This decrease of the HHG yield also results in the 
decrease of the HHG conversion efficiency [l3|, limiting 
the wider application of the HHG. 

To overcome these difficulties, great efforts have been 
devoted UMi. It has been found that the target atom 
can survive higher laser intensities in ultrashort laser 
pulses, resulting in much higher HHG cutoff energy and 
yield at high laser intensities near to the saturation ion¬ 
ization intensity (TtI-I^ . In addition, the use of two-color 
laser fields has been shown to be a very effective approach 
for increasing the HHG cutoff and producing brighter and 
shorter attosecond pulses [23 - [^ . 

The motivation of the paper is to further explore the 
procedure which could increase the HHG yield in a wide 
parameter region theoretically. We revisit the influence 
of the pulse duration on the HHG with varying the laser 
intensity and wavelength and working at both one-color 
and two-color laser fields. Unexpectedly, our simulations 
show that even for a moderate laser intensity with the low 
ionization probability, the HHG efficiency still increases 
remarkably in an ultrashort laser pulse. Our analyses re¬ 
veal a classical effect, which affects importantly on this 


phenomenon: during the fast falling part of the short 
pulse, the rescattering electron is capable of obtaining the 
same energy with traveling a shorter distance and there¬ 
fore enjoys a more efficient recollision for the HHG than it 
does in the long pulse. This classical effect becomes more 
remarkable as a two-color ultrashort pulse is used with 
increasing the HHG conversion efficiency significantly at 
diverse laser wavelengthes. Our findings have important 
implications on the dynamics of the electron in strong 
ultrashort laser pulses. 

The paper is organized as follows. In Sec. H, we in¬ 
troduce our theoretical methods. In Sec. HI, we show 
our two-dimensional (2D) numerical results for enhanced 
HHG efficiency in short laser pulses as the laser intensity 
is relatively low. The classical effect responsible for en¬ 
hancing the HHG yield is discussed in Sec. IV. Extended 
discussions for three-dimensional (3D) cases, for higher 
laser intensities and different forms of the laser envelope 
are presented in Sec. V. Sec. VI is our conclusion. 

II. THEORETICAL METHODS 

To explore a wide range of laser wavelength, we first 
use a 2D HJ model, which has been widely used in the¬ 
oretical studies [13, to simulate the HHG. 

The Hamiltonian of the model molecule studied here is 
H(t) = p^/2 -I- U(r) -b r • E(t) (in atomic units of = e = 
rrie = !)■ The potential used here has the form of V (r) = 
—Zj^^ + rl —Zj^,J^ + r 2 with 2 = (a: ± R/2cos9)‘^ + 
{y ± i?/2sind)^. R — 2 a.u. is the internuclear distance, 
^ = 0.5 is the smoothing parameter which is used to 
avoid the Coulomb singularity. Z is the effective charge 
which is modulated in such a manner that the ground 
state of the model molecule has the ionization potential of 
Ip = I.l a.u.. The latter is somewhat higher than that of 
the He atom. 9 denotes the angle between the molecular 
axis and the laser polarization. Here, we have assumed 
that the laser polarization is along the x axis and the 
molecular axis is located at the xoy plane. We consider 
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Figure 1: (Color online) Wavelength dependence of the HHG 
yield Y{N, A) for A'^-cycle one-color (a, b) and two-color (c, d) 
laser pulses. Results are obtained through full TDSE simula¬ 
tions (a, c) and the short-trajectory simulations (b, d). The 
ratio of Y[Ni)/Y{N 2 ) at A = 1400 nm is also shown. 


Figure 2: (Color online) HHG spectra for W-cycle one-color 
(a, b) and two-color (c, d) laser pulses at A = 1400 nm. Re¬ 
sults are obtained through full TDSE simulations (a, c) and 
the short-trajectory simulations (b, d), and divided by the 
cycle number N. 


the perpendicular orientation with 9 = 90*^ for which the 
molecule behaves similarly to an atom [ 2 ^. We work 
with a space grid size of x Ly = 1638.4 x 102.4 a.u. 
for the X and the y axes. The electric field used here has 
the form of E{t) = f{t)Eo{sinuJot + esm2u}ot) with e = 0 
for one-color cases and e = 0.5 for two-color cases. f{t) 
is the envelope function. To check the influence of the 
pulse duration on the HHG and simplify our discussions, 
we use a 3n-cycle laser pulse which is switched on and off 
linearly over n optical cycles with n = 1,3, 8. The whole 
pulse duration is NT with N = Zn and T = 27r/a;o. 

We solve the time-dependent Schrodinger equation 
(TDSE) using the spectral method. In each time step, a 
mask function cos^/® is used in the boundary to absorb 
the continuum wave packet. The coordinate position xq 
(?/o) from where the mask function becomes to work is 
± Lx /8 (± Ty / 8 ). Alternatively, we can set xq = ±Eq/wq 
with 2/0 = i-^y/8 unchanged [^. For our present cases, 
this treatment removes the contributions of the long tra¬ 
jectory and multiple returns to the HHG as the contri¬ 
bution of the short trajectory, which dominates experi¬ 
mental HHG [13, is not influenced basically. Eq/wq is 
the quiver amplitude of the classical electron in the laser 
field. Below, for differentiation from the full TDSE sim¬ 
ulation with xq = ±Lx/S, we denote the simulation with 
setting xq = EEq/wq as the short-trajectory simulation. 

Once the HHG power spectrum S{uj) for the harmonic 
uj is obtained from the TDSE dipole acceleration, we cal¬ 
culate the average HHG yield for a iV-cycle laser pulse 
using Y{N,X) oc /n 1 " ‘5'(w)dw. Here, Hi = Ip 


and D 2 is the cutoff energy of the spectrum. A is the 
laser wavelength. For the pulse shape used here, Y {N, A) 
can be used to compare the HHG conversion efficiency at 
different pulse durations directly M- 

In the following, our discussions will be performed for 
a moderate laser intensity of / = 5 x lO^^W/cm^ at which 
the ionization yield of the model molecule is low. 


III. ENHANCED HHG YIELD 

In Fig. 1(a), we plot the wavelength dependence of 
the average HHG yield in a one-color field for different 
pulse durations. The contrast of the curves is remark¬ 
able. One can observe that the average HHG yield is the 
highest for the 3-cycle pulse at different wavelengthes. 
This yield decreases as the cycle number N increases. At 
long wavelengthes such as A = 1400 nm, the yield of the 
3-cycle pulse is several times higher than the 9-cycle re¬ 
sult, and one order of magnitude higher than that of the 
24-cycle pulse. Note, in this case, the whole HHG yield 
of NY{N, A) for N = 3 is also several times higher than 
that for N = 24. The contrast of the HHG yields at dif¬ 
ferent pulse durations becomes more remarkable as the 
short-trajectory simulations are performed. As shown in 
Fig. 1(b), as the contributions of the long trajectory and 
multiple returns are excluded, the short-trajectory HHG 
yield for = 3 is one order of magnitude higher than 
N = 9, and two orders of magnitude higher than N = 24. 
For a two-color laser field, however, this remarkable dif- 
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Figure 3: (Color online) Rescattering time and energy dis¬ 
tributions (the color coding) for 3-cycle (a, b) and 9-cycle 
(c, d) one-color laser pulses at A = 1400 nm. Results are 
obtained through full TDSE simulations (a, c) and the short- 
trajectory simulations (b, d). In each panel, the black-square 
curve shows the electron trajectory for the first return, ob¬ 
tained from the quantum orbit model. 


our treatments and the spectrum of = 9 becomes flat, 
resulting in a remarkable contrast of these two spectra 
with different N, as shown in Fig. 2(b). 

For the two-color case, the contrast of the two spectra 
at different pulse durations is more remarkable, even for 
the full TDSE results, as shown in the right column of 
Fig. 2. In this case, the spectrum of = 3 in Fig. 2(c) 
of full simulations shows four plateaus with the cutoff 
positions of w = 1.9 a.u., 7.6 a.u., 11.3 a.u. and 17.1 
a.u., respectively. Around the second plateau which also 
has a robust peak at w = 7.6 a.u. (corresponding to 
Ep = 1.9C/p), the spectrum of Af = 3 is one order of 
magnitude higher than that of Af = 9. The latter, by 
comparison, shows two striking plateaus corresponding 
to the third and the fourth plateaus of Af = 3. As we 
perform the short-trajectory simulations, the intensity of 
the second plateau does not decrease basically, as shown 
in Fig. 2(d). The origin of the plateaus can be well under¬ 
stood through the wavelet analysis of the TDSE dipole 
acceleration combined with the quantum orbit (QO) the¬ 
ory [ 111 , 1121 , as to be discussed below. These comparisons 
in Eig. 2 explain the remarkable difference for the HHG 
yields at different pulse durations discussed in Eig. 1. 
Erom the comparisons, one can also conclude that the 
high HHG efficiency of short pulses is closely related to 
the robust HHG peaks of 1.7Up and 1.9Up indicated in 
Eig. 2, which are less influenced by our different absorb¬ 
ing procedures in TDSE simulations. This conclusion can 
be further checked through the time-frequency analysis. 


ference occurs even for the full TDSE simulations, as 
shown in Fig. 1(c). Here, the ratio of F(3, A) vs E(24, A) 
with A = 1400 nm arrives at 130, implying a significant 
increase of the HHG efficiency in a short two-color laser 
pulse. This significant increase of the efficiency is further 
amplified as the short-trajectory contributions are con¬ 
sidered, as shown in Fig. 1(d). Since the short-trajectory 
simulations are closely associated with the classical mo¬ 
tion of the electron, these results in Figs. 1(b) and 1(d) 
imply that the potential mechanism which increases the 
HHG yields at short pulses, is related to the classical as¬ 
pect of the electron. Next, we explore the mechanism in 
detail. 

Figure 2 plots the HHG spectra of 3-cycle vs 9-cycle 
pulses for the typical wavelength of A = 1400 nm. For 
comparison, these spectra are divided by the cycle num¬ 
ber N. For the one-color case in the left column of Fig. 2, 
one can observe from Fig. 2(a) of full simulations: i) the 
spectrum with A^ = 3 (the solid-black curve) is higher 
than that of A^ = 9 (the dashed-red curve), especially 
for the low-energy part, ii) In both cases, the HHG yield 
decreases as the harmonic energy increases, hi) The spec¬ 
trum of A^ = 3 shows three plateaus with the cutoff posi¬ 
tions of w = 2.5 a.u., 7.1 a.u. and 12.5 a.u. respectively. 
Around the second cutoff w = 7.1 a.u. (corresponding 
to the electron kinetic energy oi Ep = uj — Ip = 1.7t/p), 
a robust peak can be observed. As the short-trajectory 
simulation is executed, the robust peak of A^ = 3 survives 


IV. CLASSICAL EFFECTS 

As a case, in Fig. 3, we show wavelet-analysis re¬ 
sults (the color coding) for the corresponding spectra in 
Figs. 2(a) and 2(b). For comparison, here, we also show 
the electron trajectory of the first return with the excur¬ 
sion time of the electron shorter than a laser cycle (the 
black-square curve), predicted from the QO model. For 
N = 3, from Fig. 3(a) of the full TDSE results, one can 
observe that i) the distributions match the theory pre¬ 
dictions basically; ii) the distributions imply three HHG 
cutoffs. The first one with Ep = O.bUp around the return 
time tr = 2.AT has the largest amplitude. The second one 
with Ep = 1.7Up also appears around U = 2.AT and has a 
comparable amplitude with the first one, as indicated by 
the dashed arrow. The third one with Ep = 3.2Up near to 
tr = 2T has the smallest amplitude. The first one can be 
identified as arising from multiple returns, and the sec¬ 
ond and the third ones come from the first return. For 
short-trajectory simulations in Fig. 3(b), the first one 
disappears, as the second and the third ones keep their 
amplitudes with the prevailing role of the second one. 
Note, the contributions of the long trajectory around the 
minus-chirp part of the black-square curve also disappear 
basically due to the absorbing procedure used here. For 
A = 9, the situation is different. The full TDSE results 
in Eig. 3(c) show large amplitudes around the electron 
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Figure 4: The kinetic energy Ep, the maximal displacement 
Xm and the ionization time ti of the rescattering electron as 
functions of the return time tr for 3-cycle (the left column) 
and 9-cycle (the right column) one-color laser pulses at A = 
1400 nm, obtained using the quantum orbit model. Only the 
first return with the excursion time t = tr — ti < T is shown. 
The horizontal arrow indicates the orbits with Ep = 1.5Up 
for short (S) and long (LI and L2) trajectories. The vertical 
arrows indicate the corresponding Xm and ti of theses two long 
trajectories. The relevant values of Xm and ti are as shown. 


trajectories of multiple returns (these trajectories are not 
shown here), in agreement with the results in Ref. [la¬ 
in Fig. 3(d) of short-orbital simulations, the contribu¬ 
tions of multiple returns disappear. The survived dis¬ 
tributions, however, show smaller amplitudes than that 
around Ep = 1.7Up in Fig. 3(b). To understand the large 
amplitude located at Ep = l.TUp in Fig. 3(b), in Fi g. 4 , 
we further compare the maximal displacement Xm (29j, 
which the electron can travel as it ionizes at the time ti 
and returns at E tor N = 3 vs N = 9. 

To obtain the maximal displacement Xm, we first cal¬ 
culate the complex ionization time and the return time 
E (the real parts of the complex times are considered the 
physical ionization time ti and the return time tr) by the 
QO model. Then we evaluate the maximal displacement 
using dll a;™ = Xm{t) = (i?o/wQ)i?e[sma;ot - sinwo^? + 
Pst(a;g/£'o)(t-tf)] withu(t) = Re[pst+EQ/ujQCOs{u:ot)] = 
0 and Re{ti) < t < Re{t1). v{t) is the electron velocity, 
and Pst = (£;o/wo)[sina;otf -sinwot?]/[wo(i? -i?)] is the 
saddle-point momentum. We mention the maximal dis¬ 
placement obtained here agrees well with that obtained 
using the classical procedure introduced in Ref. [2^. 

Our comparisons are performed for two typical long 
trajectories (denoted using LI and L2 in Fig. 4) with 
the same return energy Ep = 1.5t/p, near to Ep = l.TUp 


N=3 N=9 
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Figure 5: Same as Fig. 4, but for two-color laser pulses. 

of the large amplitude in Fig. 3(b). Both these two 
trajectories have the ionization times ti located in the 
flat-top part of the pulse and near to the peak of the 
field, as shown in Figs. 4(c) and 4(f), and are expected 
to contribute importantly to the HHG. The return times 
of the two trajectories are different. As the Lf trajectory 
returns near to the flat-top part of the pulse, the L2 tra¬ 
jectory returns in the falling part of the pulse (the laser 
pulse becomes to fall at 2T for iV = 3 and 6T for N = 9). 
For the 3-cycle case in Fig. 4(b), the displacement Xm 
for the L2 trajectory is 97 a.u.. It is 208 a.u. for LI, two 
times larger than the L2 one. Considering the spread 
of the wave packet is proportional to the electron dis¬ 
placement, these above results imply the L2 trajectory 
has a amplitude several times larger than the LI one, in 
agreement with the wavelet-analysis result in Fig. 3(b). 
For the 9-cycle case in Fig. 4(e), the situation is differ¬ 
ent. The L2 trajectory has the maximal displacement of 
Xm = 185 a.u. which is near to 208 a.u. of the LI one. As 
a result, the wave packet spreading is comparable for the 
two trajectories in the 9-cycle case, resulting in similar 
amplitudes for them, as seen in Fig. 3(d). We mention 
that short trajectories have the smaller displacements Xm 
than the corresponding long ones. However, they usually 
ionize at a time farther away from the peak of the field 
and therefore have smaller amplitudes than the long ones 
(this can also be seen from the wavelet-analysis results 
in Fig. 3). For the reason, we don’t discuss them here. 

In combination with the distributions in Fig. 3(b) vs 
Fig. 3(d), the contrast of the maximal displacements Xm 
for the L2 trajectory in Fig. 4(b) vs Fig. 4(e) suggests 
that the higher HHG efficiency for the short pulse of A^ = 
3, observed in Fig. 1(b), is closely related to the shorter 
excursion distance of the electron in the fast falling part 
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of the short pulse. This classical effect which increases 
the HHG efficiency becomes more remarkable in the two- 
color case, as shown in Fig. 5. Here, our analyses are also 
performed for two typical long trajectories of LI and L2. 
The ionization times of the two long trajectories both 
are located at the flat-top part of the pulse and are near 
to the peaks of the two-color field, as shown in Figs. 
5(c) and 5(f). In addition, the electric-field amplitudes 
at the two ionization times of LI and L2 in the two- 
color case are nearer to each other than in the one-color 
one. For the short-pulse case of = 3, one can observe 
from Fig. 5(b) that the maximal displacement of the L2 
trajectory is Xm = 77 a.u., as the LI trajectory shows 
a maximal displacement of Xm = 217 a.u., almost three 
times larger than the L2 one. For the long-pulse case of 
N = 9, they are 154 a.u. and 217 a.u., respectively, as 
shown in Fig. 5(e). The shorter excursion distance of 
the L2 trajectory in Fig. 5(b) of fV = 3, suggests that 
this trajectory contributes significantly to the HHG in 
the short-pulse case. It is corresponding to the cutoff 
of the second plateau with high intensity in Figs. 2(c) 
and 2(d). We mention that the TDSE cutoff position 
of the second plateau in Fig. 2(c) is 1.9Up for the 3- 
cycle case, somewhat higher than the model prediction 
of 1.7Up in Fig. 5(a). This is also the one-color case 
of 1.7Up in Fig. 2(a) versus l.bUp in Fig. 4(a). This 
difference can partly arise from the nonadiabatic effect in 
ultrashort pulses. In comparison with Xm = 97 a.u. with 
Ep = 1.5Up in Fig. 4(b), this shorter excursion distance 
Xm = 77 a.u. of the L2 trajectory with the higher return 
energy Ep = 1.71/p in Fig. 5(b) also suggests the HHG 
efficiency is higher in the two-color case than in the one- 
color case, in agreement with our analyses in Fig. 1. 

V. EXTENDED CONSIDERATIONS 
A. Three-dimensional simulations 

To check our results, we have also performed 3D sim¬ 
ulations for HJ with the soft-core potential of V{r) = 
+r\ — -I- r| and 2 = (a: ± i?/2cos0)^ -I- 

{y ± i?/2sin0)^ -|- . The Definitions of the parameters 

Z, R and 9 are the same as in our 2D cases. The 3D 
calculations are very time-memory consuming. Here, we 
work with a grid size of L^x LyX L^ = 819.2 x 51.2 x 51.2 
a.u. for the x, y and z axes, respectively. Our calcula¬ 
tions are performed for / = 5 x 10^'^W/cm^ and A = 1400 
nm corresponding to the laser parameters used in Fig. 2. 
Similar absorbing procedures as in 2D cases are also used 
in performing full TDSE simulations and short-trajectory 
simulations. The results are presented in Fig. 6. 

One can observe from Fig. 6 that the 3D results are 
similar to the 2D ones in Fig. 2. First, the spectra of 
N = 3 show a robust peak which appears around Ep = 
oj — Ip = 1.7Up for one-color cases in the left column of 
Fig. 6 and around Ep = 1.91/p for two-color cases in 
the right column of Fig. 6. Secondly, this robust peak 
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Figure 6: (Color online) Same as Fig. 2, but obtained with 
3D simulations. 

is more remarkable in short-trajectory simulations than 
in full simulations and in two-color cases than in one- 
color cases. As discussed in Fig. 2, this robust peak 
arises from the classical effect, which increases the HHG 
efficiency importantly in short pulses. The 3D results in 
Fig. 6 also show that this HHG efficiency is strikingly 
higher in the short pulse than in the lone one. All of the 
characteristics are in agreement with our 2D results in 
Fig. 2. 

One of the main differences between 2D and 3D cases is 
that the diffusion effect relating the wave packet spread¬ 
ing is stronger in 3D cases than in 2D cases. This stronger 
diffusion effect can induce a larger phase difference be¬ 
tween the harmonics emitted in different laser cycles, es¬ 
pecially for harmonics arising from the long trajectory 
and multiple returns with longer excursion times in the 
laser field. The interference of the harmonics emitted 
in different laser cycles will decrease the intensity of the 
HHG spectra and this decrease is more remarkable for 
long pluses with more cycles. This can be the reason that 
the spectrum of A = 9 in Fig. 6(a) of full TDSE simula¬ 
tions shows a smaller amplitude in comparison with that 
of A = 3 in the high-energy region with ui > 7 a.u.. Note 
that in Fig. 6(b) of short-trajectory simulations, these 
two spectra of A = 3 and A = 9 are comparable for the 
high-energy region of w > 7 a.u., as the contributions of 
the long trajectory and multiple returns are excluded. 

B. Sin-square-envelope pulses 

To check our results, a sin-square-envelope pulse is also 
used in our calculations and relevant results are presented 
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Figure 7: (Color online) Same as Fig. 2, but obtained for a 
sin-square-envelope laser pulse. 

in Fig. 7. The results in Fig. 7 are also similar to those 
in Fig. 2, with showing a robust peak in the spectra of 
N = 3. Here, the position of the peak is around Ep = 
O.SUp for one-color cases and Ep = Up for two-color cases, 
somewhat lower than those in Fig. 2. In addition, the 
maximal cutoff position of the spectrum of = 3 is 
also somewhat lower than that of = 9. However, in 
comparison with the results in Fig. 2, the results for 
N = 3 presented here show higher inversion efficiency. 
For example, in Fig. 7(b), the spectrum of A^ = 3 is 
higher than that of A^ = 9 in the whole energy region, 
different from the results in Fig. 2(b). In addition, the 
curve of A^ = 3 in Fig. 7(d) also shows a cutoff located at 
w = 13.5 a.u.. Around the cutoff, the spectrum of A^ = 3 
is two orders of magnitude higher than that of Af = 9. 
All of the characteristics can be understood by virtue 
of the analyses of the quantum orbit and the maximal 
displacement of the rescattering electron. 

C. High laser intensities 

As the laser intensity increases and is near to the sat¬ 
uration intensity, the situation is somewhat different. As 
shown in Fig. 8(a), for the case of / = 1.2 x 10^®W/cm^ 
and A = 1400 nm, the HHG spectrum of A^ = 3 is higher 
than that of A^ = 9 in the whole energy region, even the 
average over the cycle number N is not performed here. 
This can be understood from the ground-state deletion. 
For the long pulse of A^ = 9, the ionization is stronger 
than that of A^ = 3, resulting in a significant deletion of 
the ground state and accordingly a remarkable decrease 
of the HHG yield. However, in this case, the classical 
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Figure 8: (Color online) HHG spectra for A-cycle one-color 
laser pulses at 7 = 1.2 x 10^®W/cm^ and A = 1400 nm (a, 
b) and 7 = 1 x 10^®W/cm^ and A = 1600 nm (c, d). Results 
are obtained through full TDSE simulations (a, c) and the 
short-trajectory simulations (b, d) in 2D cases. 

effect can still be read from the spectrum, which mani¬ 
fests itself as a harmonic peak around Ep = 1.75t/p, as 
indicated by the dashed-dotted arrow. This peak is more 
remarkable in short-trajectory simulations, as shown in 
Fig. 8(b). In addition, the cutoff position in the spec¬ 
trum of A^ = 3 in Fig. 8(a) seems somewhat larger 
than that of A^ = 9. This phenomenon is clearer for 
the case of 7 = 1 x lO^^W/cm^ and A = 1600 nm, as 
shown in Fig. 8(c). We mention that for the case in 
Fig. 8(c), the ionization is somewhat weaker than that 
in Fig. 8(a). Accordingly, the spectra of Af = 3 and 
N = 9 are comparable here. In addition, a harmonic 
peak around uj = 1.787/p can also be observed in the 
spectrum of A^ = 3 in Fig. 8(c). This peak is more strik¬ 
ing in Fig. 8(d) of short-trajectory simulations. These 
results show that the classical effect, mainly discussed in 
the paper, is general in a wide laser-parameter region. 


VI. CONCLUSIONS 

In summary, we have shown that the HHG efhciency 
in an ultrashort laser pulse is influenced significantly by 
a classical effect. The latter is closely associated with the 
shorter excursion distance of the rescattering electron as 
it ionizes near the peak of the short laser pulse and re¬ 
turns in the fast falling part of the pulse. This shorter ex¬ 
cursion distance suppresses the spread of the wave packet 
and increases the efficiency of the HHG. With the shorter 
excursion distance and accordingly the shorter excursion 
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time, the emitted harmonics relating to the classical ef¬ 
fect can be easier to survive the macroscopic propagation 
in the medium and easier to modulate in experiments. 
We expect that this effect also has an important influence 
on the rescattering induced other strong-held phenomena 
such as high-order above threshold ionization and nonse¬ 
quential double ionization. 
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